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ABSTRACT 

The HH 111 protostellar system is a young Class I system with two sources, 
VLA 1 and VLA 2, at a distance of 400 pc. Previously, a flattened envelope has 
been seen in C 18 to be in transition to a rotationally supported disk near the 
VLA 1 source. The follow-up study here is to confirm the rotationally supported 
disk at 2—3 times higher angular resolutions, at ~ 0!'3 (or 120 AU) in 1.33 mm 
continuum, and ~ 0"6 (or 240 AU) in 13 CO (J = 2 - 1) and 12 CO (J = 2 - 1) 
emission obtained with the Submillimeter Array. The 1.33 mm continuum emis- 
sion shows a resolved dusty disk associated with the VLA 1 source perpendicular 
to the jet axis, with a Gaussian deconvolved size of ~ 240 AU. The 13 CO and 
12 CO emissions toward the dusty disk show a Keplerian rotation, indicating that 
the dusty disk is rotationally supported. The density and temperature distribu- 
tions in the disk derived from a simple disk model are found to be similar to those 
found in bright T-Tauri disks, suggesting that the disk can evolve into a T-Tauri 
disk in the late stage of star formation. In addition, a hint of a low-velocity 
molecular outflow is also seen in 13 CO and 12 CO coming out from the disk. 

Subject headings: circumstellar matter - stars: formation — ISM: individual (HH 
111) 

1. Introduction 

Stars are formed inside molecular cloud cores by means of gravitational collapse. The 
details of the process, however, are complicated by the presence of magnetic fields and angular 
momentum. In particular, for the stars to form, collapsing material has to lose most of the 
magnetic fields and angular momentum. As a result, in addition to infall (or collapse), 
rotation and outflow are also seen toward star-forming regions. In theory, a rotationally 
supported disk is expected to form inside a collapsing core around a protostar, from which 
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part of the material is accreted by the protostar and part is ejected away. Observationally, 
however, when and how a rotationally supported disk is actually formed are still unclear, 
because of the lack of detailed kinematic studies inside the collapsing core. 

This paper is a follow-up study of the HH 111 protostellar system, which is a well-studied 
Class I system deeply embedded in a compact molecular cloud core in the L1617 cloud of 
Orion at a distance of 400 pc. This system has a bolometric temperature T hn i ~ 7 8 K and 



thus has just transitioned to the Class I phase from the Class phase (jFroebrichl 120051 ). 
At the center of this system, t here are two sources, VLA 1 and VLA 2, with a projected 
separation of ~ 3" (1200 AU) (IReipurth et al.lll999l ). Previously, a flattened envelope has 
been seen in C 18 around the VLA 1 source (Lee et al 2009, hereafter Paper I) and it seems 
to be in transition to a rotationally supported disk near the VLA 1 source (Lee et al 2010, 
hereafter Paper II). The follow-up study here is to confirm the rotationally supported disk 
at 2—3 times higher angular resolutions, at ~ 0''3 (or 120 AU) in 1.33 mm continuum, and 
~ 0"6 (or 240 AU) in 13 CO (J = 2-1) and 1 2 CO (J = 2 - 1) emission obtained with 
the Submillimeter Array (SMA)0 (jlio et al. 2004 ). This study also shows that the physical 
properties of the disk here in the early phase of star form ation can be similar to those of the 
bright T-Tauri disks in the late phase of star formation ([Andrews et al.l 120091 ). 



2. Observations 

Observations toward the HH 111 system we re carried out with the SMA, first on 2005 



December 5 in the compact-north configuration (IPaper J ), then on 2010 Januar y 24 in th e 
extended configuration and on 2010 March 26 in the sub-compact configuration (IPaper HI ), 
and then on 2011 Jan 01 in the very extended configuration. The baselines have projected 
lengths ranging from ~ 10 to 480 m onto the sky. The 230 GHz band receivers were used to 
observe the 12 CO (J = 2-1), 13 CO (J = 2-1), C 18 (J = 2-1), and SO (Nj = 5 6 -4 5 ) lines 
simultaneously with the 1.33 mm continuum. The 12 CO and 13 CO lines are bright and thus 
are presented here to show the kinematics of the inner region and search for the rotationally 
supported disk around the VLA 1 source at high resolution. The velocity res olution in these 



lines is ~ 0.28 km s -1 per channel, similar to that of C 18 presented in IPaper IB . The 
visibility data were calibrated with the MIR package. The flux uncertainty was estimated 
to be ~ 15%. The calibrated visibility data in different configurations were combined and 
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then imaged with the MIRIAD package, as described in iPaper II With various different 
weightings, the synthesized beams can have a size from ~ O'.'Q down to 0"3. The rms noise 
levels are ~ 0.036 Jy beam -1 in the 12 CO channel maps with a beam of 0"58x0"50, 0.032 Jy 
beam -1 in the 13 CO channel maps with a beam of 0!'66x0"56, and 0.9 mJy beam -1 in the 
continuum map with a beam of 0''36x0''31. The velocities in the channel maps are LSR. 



Results 



As in IPaper III , the results are presented in comparison to a mosaic image based on the 
Hubble Space Tel escope (HST) N I CMO S image ([Fell] 1.64 fim + H 2 at 2.12 fim + contin- 
uum) obtained by lReipurth et al.l (11999 ). which shows the two sources VLA 1 and VLA 2 in 
the infrared, and the reflection nebulae that trace the illuminated outflow c avity walls. The 



two sources are located at a 



2008) and 



a 



(20001 



(2000) 

05 h 51 m 46?07, 5 



05 h 51 m 46?254, 6, 



(2000) 



(2000) 



+02°48 / 30: / 76, respectively flReipurth et al 



+02°48 / 29 "65 ([Rodriguez et al. 



1999 



Rodriguez et al.ll2008l ). The jet is bipolar with a western and an eastern component. It has a 
position angle (P. A.) of ~ 97° and an inc lination angle of ~ 10°, with the wes tern component 



being blueshifted and tilted toward us (IReipurth. Raga. &: Heathcotdll992l ). The systemic 
velocity in this system is assumed to be 8.90±0.14 km s -1 LSR, as before. Throughout this 
paper, the velocity is relative to this systemic value. 



3.1. 1.33 mm Continuum Emission 



As shown in Figure [T^, the continuum emission associated with the VLA 1 source is now 
resolved into a disk-like structure perpendicular to the jet axis, with a Gaussia n deconvo lved 
(FWHM) size of ~ 0"6x0"3 (~ 240x120 AU) and a P.A. of 7.6°. As found in lPaper ij this 
emission has a total flux of ~ 285±40 mJy and it is the thermal emission mainly from a 
dusty disk, wit h the inner part of whi ch already seen with a deconvolved size of ~ 0''15 (60 
AU) at 7 mm ([Rodriguez et al.l 120081 ) . Note that the emission seems to show a little faint 
protrusion to the east along the jet axis. This faint protrusion seems to be seen in 3.6 cm 
and 7 mm as well, and it may trace either the mater ial along the jet axis or a disk around a 
companion in a close binary system as suggested in iRodrfguez et al.l (120081 ) . In addition, a 
faint protrusion is also seen extending to the southeast from the southern part of the disk. 
As in IPaper III , faint emission is also detected around the VLA 2 source (~ 4 a detection 
with lcr = 0.9 mJy beam -1 ), probably tracing a disk around it. Moreover, two more faint 
emission peaks DCN and DCS are seen at 3" in the north and 2" in the south respectively 
near the equatorial plane (P.A. ~ 7°), probably tracing the density enhancement there. 
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3.2. Previous Results of Envelope and Outflow in C ls O 



Previously in the C O observations ([Paper It ), a flattened envelope is seen perpen- 



dicular to the jet axis, extending to ~ 7200 AU (18") out from the VLA 1 source. It is 
rotating, with the blueshifted emission in the north and the redshifted emission in the south. 
The outer part (~2000— 7200 AU, or ~ 5"-18") seems to have a rotation that has constant 
specific angular momentum with ~ 3.90<i _1 km s , while the inner part (~ 60—2000 AU, 
or ~ 0''15-5") seems to have a Keplerian rotation with v$ ~ 1.75c? -0 " 5 km s _1 , where ci is the 
radial distance from the source in the unit of arcsecond. Outflow is also seen in C 18 0, with 
the blueshifted emission extending to the west and the redshifted emission extending to the 
east. 



3.3. Line Emission 

The 13 CO and 12 CO lines are detected toward the VLA 1 source upto ~ 4.2 km s -1 and 
10 km s _1 , respectively, from the systemic velocity. In order to show how the envelope and 
disk structure changes with velocity, the line emission is divided into four velocity ranges: low 
(0—2.5 km s _1 ), medium (2.5—4.2 km s _1 ), high (4.2—5.8 km s -1 ), and very high (5.8—10 
km s _1 ) velocity ranges, on the redshifted and blueshifted sides. 

In 13 CO, at low velocity, the redshifted emission is detected mainly in the south and 
the blueshifted emission is mainly in the north of the VLA 1 source (Fig. [2h), consistent 
with a rotation motion about the source. As seen in C 18 0, the emission is contaminated 
by the outflow emission, with the blueshifted emission also extending to the west and the 
redshifted emission also extending to the east. In addition, the blueshifted emission shows 
a V-shaped structure opening to the north, spatially coincident with the cavity walls traced 
by the reflection nebulae, suggesting that the envelope material is piling up on the cavity 
walls as the cavity walls expand into the envelope. A similar V-shaped structure is also seen 
opening to the south in the redshifted emission. At medium velocity, the 13 CO emission 
shrinks to the source, spatially coincident with the dusty disk (Fig. [2^. and Fig. [Tb for a 
zoom-in). The blueshifted emission peak is to the north and the redshifted emission peak is 
to the south in the equatorial plane, indicating that the motion there is highly dominated by 
the rotation about the VLA 1 source, as expected for a rotationally supported disk. The two 
peaks, however, are not symmetric about the source, with the northern peak at a distance 
of ~ 0"25 and southern peak at ~ 0''5. Two faint protrusions are also seen, one in the blue 
extending to the west from the northern part of the disk and one in the red extending to the 
east from the southern part, parallel to the jet axis. These structures are similar to those 



seen in the rotating molecular outflow in CB 26 (ILaunhardt et al.ll2009l ). and thus may trace 
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material outflowing from the disk. 

In 12 CO, at low velocity, the blueshifted emission is mainly to the north and the red- 
shifted emission is mainly to the south, with V-shaped structures coincident with the cavity 
walls (Fig. [2b), similar to that seen in 13 CO. At medium velocity, the blueshifted emission 
extends to the west from the northern part of the disk while the redshifted emission extends 
to both the west and east from the southern part of the disk (Fig. [2H and Fig. [lb for a zoom- 
in). Note that the blueshifted emission and the redshifted emission beyond ~ ±1" from the 
source are contaminated by the outflow shell emission, showing a cone-like structure with 
the tip pointing toward the source. The blueshifted emission also shows an E-W elongation 
at ~ 3" in the north (Fig. [2bl), with an unknown origin. At high velocity, the blueshifted 
and redshifted emissions extend out from the inner radii of the disk (Fig. [2b and Fig. [Tbl 
for a zoom- in). The emissions within ~ 1" from the source are not much contaminated by 
the outflow shell emission, which is expected to be further away from the source at higher 
velocity. They extend to the west, first opening from the VLA 1 source and then bent to 
be aligned with the jet axis. They also extend slightly to the east. The blueshifted emission 
is mainly to the north of the jet axis and the redshifted emission is mainly to the south, 
with the velocity sense the same as that of the rotation in the disk, and thus may trace the 
material coming out from the disk. Note that, in the west, the blueshifted emission also 
extends slightly to the south. At very high velocity, only blueshifted emission is detected 
with a peak slightly away from the equatorial plane at ~ 0''25 to the west (Fig. [2]F and Fig. 
[Tb for a zoom- in). This emission peak is spatially coincident with the blueshifted emission 
peak seen at high velocity in Figure [Tbl, suggesting that the emission there has a broad range 
of velocities from high to very high and may trace a shock interaction, either inside the disk 
or outside. 



3.3.1. Rotation motion near the VLA 1 source 

As seen in the last section, the inner part of the envelope can also be traced by 13 CO 
and 12 CO, allowing us to check and confirm the Keplerian rotation law near the source seen 
in C 18 0. Figure [3] shows the position- velocity (PV) diagrams of the envelope and disk in 
13 CO and 12 CO cut along the equatorial plane perpendicular to the jet axis. In 13 CO, the 
emission within ±1 km s -1 from the systemic velocity is mainly resolved out and affected 
by the outflow emission, and thus should be excluded for studying the rotation law. As can 
been seen, the velocity structure beyond ±1 km s _1 roughly follows the Keplerian rotation 
law seen in C 18 0. In 12 CO, the emission within ±2 km s _1 from the systemic velocity is 
mainly resolved out and affected by the outflow emission, and thus should also be excluded 
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for studying the rotation law. In addition, as seen in the emission map (Fig. |2ji and e), the 
blueshifted emission is also detected at ~ 2" to 3" in the north with a velocity ranging from 
—2 to —6 km s _1 w.r.t. the systemic velocity. This emission is unclear, and it may trace 
some interaction in the envelope or something in the foreground or background, and thus 
should also be excluded. In this line, on the redshifted side, the emission is detected up to 
~ 6 km s _1 , higher than that seen in 13 CO, arising from the inner region as seen in Figure 
[T]i. The velocity structure also follows the Keplerian rotation law to that high velocity to 
the inner region. On the blueshifted side, the emission is detected up to —10 km s -1 from 
the systemic, much higher than that on the redshifted side. However, as can be seen, the 
emission beyond —5 km s _1 is localized with a broad range of velocities and it corresponds 
to the localized emission at ~ 0''25 to the west as seen in the emission map (Fig. [Tfe). Since 
this emission can trace a local shock interaction, the blueshifted emission should not be used 
to derive the rotation law. As a result, the rotation structure near the source seen in 12 CO 
is also consistent with that seen in C 18 0. 



Model the disk continuum emission 



Here a simple disk model is used to fit the structure and intensity distribution of the 
continuum emission in order to derive the temperature and density distributions of the dusty 
disk. Note that, the emission that extends to the east along the jet axis and the emission 
that extends to the southeast from the southern part of the disk are not from the disk 
perpendicular to the jet axis and are thus excluded from the model. 

For simplicity, the disk in the model is assumed to be flat with a constant thickness of H, 
an inner radius of i? in , and an outer radius of R ou t- The disk is assumed to be perpendicular 
to the jet axis an d thus have an inclination angle of 1 0° to the line of sight, with the nearside 
tilted to the east (IReipurth. Raga. &: Heathcotdll992l ). It has a position angle of 7.6° as given 
by the Gaussian source deconvolution earlier in Section 13.11 The temperature is assumed to 
be given by 

nn / R 



T 



Rn 



1 



where Rq = 0''3 (120 AU), which is the radius (i.e., the half of the size) of the disk derived 
earlier from the Gaussian source deconvolution, T is the temperature at Rq, and the power- 
law index q is set to 0. 5, assumed to be s i milar to that found in T -Tauri d isks in the late 
stage of star formation (IMundy et al.lll996t Andrews et al.ll2009l ). In lPaper J . using the dust 
opacity law 

0.1 f ,,„^„ Y cm 2 g" 1 (2) 



10 12 Hz 
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as given in lBeckwith et al.l (119901 ) . the spectral energy distribution (SED) of the continuum 
source was found to be fitted with a dust temperature ranging from 41 to 64 K and /3 ~ 0.9. 
Thus, here T is assumed to be 45 K and the dust mass opacity is assumed to be 0.026 cm 2 
g _1 at A = 1.33 mm (or 225 GHz). The number density of molecular hydrogen is assumed 
to be given by 

n = no(^r)- p (3) 

with the power-law index p assumed to be 1. With a constant thickness, the disk has a 
surface density with a power- law index of 1, similar to that found in the T-Tauri disks 
( [Andrews et al.l 120091 ) . As can be seen in Figure HI with the above assumptions, the model 
can fit the continuum emission reasonably well, with H ~ 0''3 (120 AU), R[ n ~ 0!'03 (12 AU), 
-Rout ~ 0!'6 (240 AU), and n ~ 10 9 cm -3 . Here, H turns out to be the same as the Gaussian 
source deconvolved size of the continuum emission in the minor axis. This is expected for 
a nearly edge-on disk, for which the size in the minor axis is mainly due to the thickness of 
the disk. On the other hand, R out ~ 2Rq and thus is twice the disk radius derived from the 
Gaussian source deconvolution. Assuming that Helium has a number density % e = 0.2n, 
then the mean molecular weight with respect to the number of molecular hydrogen is 2.8. 
Thus, the mass of the di sk is found to be 0.14 M , similar to the mean value found by 
fitting the SED (IPaper it ). The surface density in the disk is ~ 7.2 g cm -2 at 120 AU and 
21.5 g cm -2 at 40 AU, similar to those found in the bright T-Tauri disks, e.g., GSS 39, in 
Ophiuchus, which als o has a similar radius of ~ 200 AU and a same disk mass of ~ 0.14 M 
([Andrews et al.l 120091 ) . Thus, the dusty disk here can be a very young protoplanetary disk 
that could evolve into a T-Tauri disk in the later stage of star formation. 



5. Discussion 
5.1. A Rotationally supported Disk 

A rotationally supported disk has been seen in Taurus with an outer radius of ~ 500 AU 



in CO in the late (i.e., Class II or T Tauri) phase of star for mation (see, e 



Simon et al. 



20001 ). Note that the disk could appear smaller in continuum ([Andrews et al.ll2009l ). because 
of the rapid decrease in the continuum intensity due to a power-law decrease of the density 
and temperature in the disk. Formation of such a disk must have started early in the Class 
I phase and even earlier in the Class phase. HH 111, with a bolometric temperature of 
only ~ 78 K, is a very young Class I source, just transitioning from the Class phase. 
The observations here show that such a rotationally supported disk can indeed form in such 
an early phase of star formation. In continuum, the disk is found to have a deconvolved 
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radius of ~ 120 AU (with an outer radius of 240 AU in the disk model) , with the inner 
part of which already seen with a deconvolved radius of ~ 30 AU at 7 mm ([Rodriguez et al. 



20081 ). It is rotationally supported with a Keplerian rotation, as seen previously in C 18 an d 
now in 13 CO and 12 CO, which im plies a mass of ~ 1.3 M for the central star (IPaper it ). 
Note that as discussed in IPaper III , this stellar mass can account reasonably for the observed 
bolometric luminosity. The disk mass is ~ 0.14 M , which is ~ 10% of the stellar mass. This 



disk mass is consistent with that implied for a Class source by Jjprgensen et al.l ( 120091 ) and 
that implied for a young Class I source by Andrews fc Williams! ( 120071 ). Since no rotationally 
supported disk has been confirmed in the Class phase, the disk here could be the youngest 
rotating disk resolved to date. Moreover, as discussed earlier, the surface density and the 
radius of the disk here in HH 111 a re both similar to those of the bright resolved T-Tauri 
dusty disks in Andrews et al.l (120091 ) . Thus, the disk could evolve to a T-Tauri disk in the 
later stage of s tar forma tion as the accretion rate drops to 10~ 7 M yr _1 from a few times 
10~ 6 M yr" 1 JPaper III ). 



As discussed in IPaper it the Keplerian rotation in this system can extend further out 
to the envelope to ~ 2000 AU (5") from the source. At that large radius, the envelope 
seems to have some small amount of infall motion as well, and thus the rotation there could 
actually be sub-Keplerian. Therefore, the disk seems to be deeply embedded in a sub- 
Keplerian envelope with a radius of ~ 2000 AU. Observations at higher velocity and angular 
resolutions are really needed to differentiate the Keplerian motion from the sub-Keplerian 
motion in order to determine where the envelope is transitioning to a disk. Another way 
to determine the location of transition is t o search for an accretion shock that is expected 
to occur at the transition. As discussed in IPaper III . SO is found to trace a shock emission 
at 400 AU (1")) and it may trace an accretion shock and thus the transition there. Further 
observations in other shock tracers are needed to check this possibility. 

Rotationally supported disks have also been claimed in a few other Cla ss I sources, 
e.g., IRAS 04302 +2247 (the Butter fly star, T bol = 203 K) rtwblf et al.l 120081) L1489-IRS 
(T bol = 238 K) (IBrinch et al.ll2007h. IR S 43 (T bol = 310 K) JjQrgensen et al.l I2OO9I ) . and 
Elias 29 (T bo i = 350 K) (ILommen et al.ll2008l ). These sources, with a much higher T bo i than 
HH 111, are at much later stages than HH 111. Among these sources, L1489-IRS has been 
modeled in detail for the disk and envelope properties, and thus can be compared here with 
HH 111. This system was found to host a Keplerian disk with an outer radius of ~ 200 AU 
deeply embedded in a sub-Keplerian envelope with a radius ~ 2000 AU (IBrinch et al.ll2007l ). 



similar to the HH 111 system. The mass and the inf all rate in the sub-K eplerian envelope are 
~ 0.093 M & and ~ 4.3xl0~ 6 M yr _1 , respectively ( Brinch et al. 20071), also si milar to those 
in HH 111, which are ~ 0.13 M and 4.2xl0 -6 M yr" 1 , respectively ( IPaper it ). The central 
star has a mass of ~ 1.35 M , also similar to that in HH 111. However, this system has a 
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bolometric luminosity of ~ 3.7 L , much smaller than that of HH 111, which is ~ 20 L & . The 
small bolometric luminosity suggests that the accretion rate in the disk must be much lower 
than that in HH 111, as supported by the less jet activity. This is probably because that the 



only 4x10" 



.8 Mr, ( 


Joreensen et al. 


Jorgensen et al. 


(2009 





however, it is not clear if this is because the disk mass really decreases with the evolutionary 
stage, or merely because the disk mass is underestimated more in the later stage due to a 
grain growth. 

Rotationally supported disks have b een seen in a few other much older s ources, eg., HH 



30 (late Class I) (IGuilloteau et all 120081 ) and IRAS 04158+2805 (Class II) (Andrews et al. 
20081 ). In these older disks, a central hole is seen and it could be due to a tidal clearing by a 
possible binary system at the center or due to an outflow cavity. Since the HH 111 VLA 1 
source itself is a close binary, a central hole is expected to have formed at the center of the 
disk near the source as well. Observations at higher angular resolution are needed to check 
this possibility. 



5.2. Material outflowing from the disk? 

As mentioned in Section 13. 3[ 12 CO and 13 CO emissions are seen near the source ex- 
tending out from the disk to the east and west parallel to the jet axis. Co uld this emission 



struc ture trace the material coming out from the disk, as that seen in CB 26 (ILaunhardt et al. 



20091 ). carrying away the angular momentum? To study this possibility, PV diagrams cut 
across the jet axis centered at 5 positions along the jet axis, from — 1'/2 in the east to 1('2 in 
the west with a step of 0''6 (which is the size of the synthesized beam), are presented here 
in Figure [5] to investigate the kinematics of this emission structure. In these cuts, excluding 
the emission within ±2 km s -1 from the systemic velocity, the blueshifted emission is in 
the north and redshifted in the south, as seen in the emission map in Figure [TU, with a 
velocity sense the same as that of the rotation in the disk. Therefore, the emission could 
trace the material rotating around the jet axis and thus the material outflowing from the 
disk. In panels (c) and (d), the blueshifted emission is contaminated by the shock emission 
and should be excluded. From cut A in the east to cut E in the west, no clear systemic 
velocity shift is seen in the PV diagrams. Thus, if the emission really traces the material 
outflowing from the disk, the projected outflow velocity must be smaller than 1 km s" 1 . This 
is not impossible since (1) the outflowing material could have the same inclination as the 
jet and thus be close to the plane of the sky with an inclination angle of ~ 10°, and (2) the 



outflowing velocity could be only a few km s 1 in 12 C0 ( ILaunhardt et al 
observations are needed to study this possibility. 



20091 ). Further 



6. Conclusion 

The HH 111 protostellar system has been mapped at a resolution of ~ 0!'3 in 1.33 mm 
continuum, and 0!'6 in 13 CO (J = 2 — 1) and 12 CO (J = 2 — 1). The 1.33 mm continuum 
emission shows a resolved dusty disk associated with the VLA 1 source perpendicular to 
the jet axis, with a Gaussian deconvolved size of ~ 240 AU. The 13 CO and 12 CO emissions 
toward the dusty disk show a Keplerian rotation, indicating that the dusty disk is rotationally 
supported. The density and temperature distributions in the disk derived from a simple disk 
model are found to be similar to those found in bright T-Tauri disks. Thus, the disk can be 
a young protoplanetary disk that could evolve into a T-Tauri disk in the late stage of star 
formation. In addition, a hint of a low- velocity molecular outflow is also seen in 13 CO and 
12 CO coming out from the disk. 

I thank the SMA staff for their efforts in running and maintaining the array. I also 
thank the referee for the valuable comments. 
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Fig. 1. — 1.33 mm continuum, 13 CO, and 12 CO emiss ion contours overplotted on the HST 
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NICMOS image adopted from iReipurth et al.l (Il999[ ). The asterisk and cross mark the 
positions of the VLA 1 and 2 sources, respectively. The blue and red arrows indicate the 
orientations of the blueshifted and redshifted components of the jet, respectively. The dashed 
line indicates the equatorial plane perpendicular to the jet axis. The synthesized beam has 
a size of 0"36x0"31 in continuum, 0"66x0"56 in 13 CO, and 0"58x0!'50 in 12 CO. The gray 
contours are the continuum contours with levels of 4<r(l — r n )/(l — r), where r = 1.45, 
n=l,2,3.., and a = 0.9 mJy beam -1 , (b) shows the redshifted (red contours) and blueshifted 
(blue contours) 13 CO emission at medium velocity. The contour levels start at 4cr with a 
step of 3a, where a = 0.023 Jy beam -1 , (c), (d), and (e) show the redshifted (red contours) 
and blueshifted (blue contours) 12 CO emission at medium velocity, high velocity, and very 
high velocity, respectively. The contour levels start at 4cr with a step of 6a in (c), and a step 
of 4<T in (d) and (e), where a = 0.025 Jy beam -1 . 
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Fig. 2. — Same as Figure [T] but adding low-velocity 13 CO and 12 CO emission contours, 
without zooming into the central region. Here (b) corresponds to Figure [lb, and (d)-(f) 
correspond to Figure Hb-e. (a) shows the the redshifted (red contours) and blueshifted (blue 
contours) 13 CO emission at low velocity. The contour levels start at 4cr with a step of 6a, 
where a = 0.023 Jy beam" 1 , (c) shows the the redshifted (red contours) and blueshifted 
(blue contours) 12 CO emission at low velocity. The contour levels start at 4<r with a step of 
8cr, where a = 0.025 Jy beam -1 . 
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Fig. 3. — Position-velocity (PV) diagrams in (a) 13 CO and (b) 12 CO centered at the VLA 
1 source cut along the equatorial plane perpendicular to the jet axis. The gray box in the 
lower-right corner shows the velocity and spatial resolutions of the PV diagrams. Solid curve 
shows the rotation that has constant specific angular mome ntum and dashed curve shows 
the Keplerian rotation, as derived from the C 18 envelope in |Paper~Til . with both corrected 
for the inclination angle. Here d is the radial distance from the source in arcsec. In (a), the 
contour levels start at 3a with a step of 3a, where a = 0.023 Jy beam -1 . In (b), the contour 
levels start at 3a with a step of 4a, where a = 0.025 Jy beam -1 . 



- 15 - 




Fig. 4. — A fit to the continuum emission with a simple disk model as described in Section 
HI The asterisk marks the position of the VLA 1 source. The contour levels are the same as 
those in Figure [Jji. 
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Fig. 5. — PV diagrams in 13 C0 (left column) and 12 C0 (right column) cut across the jet 
axis centered at 5 positions along the jet axis, from (a) —1"2 in the east to (e) 1"2 in the 
west with a step of 0!'6 (which is the size of the synthesized beam). The gray box in the 
lower-right corner of (a) shows the velocity and spatial resolutions of the PV diagrams. The 
solid and dashed curves, and contour levels are the same as those in Figure [3j 



